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Possible causes of a rise with energy of the cosmic ray positron fraction 
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Based on general considerations rather than model-dependent fits to specific scenarios, we argue 
that an increase with energy of the positron fraction in cosmic rays, suggested by several experiments 
s-t -E ^ 7GeV, most likely requires a primary source of electron-positron pairs. We discuss the 
possible alternatives, and find none of them plausible on astrophysical or particle physics grounds. 
Further observational ways to test different scenarios are discussed. 



PACS numbers: 98.70.Sa 
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I. INTRODUCTION 

Since longtime, the study of the positron/electron ratio 
in cosmic rays has been recognized as an important tool 
to constrain the production and propagation of energetic 
particles in the Galaxy and in the Solar System. The 
PAMELA satellite detector, which began its three-year 
mission in June of 2006, is designed to measure (among 
other components) the spectra of cosmic ray positrons up 
to 270 GeV and electrons up to 2 TeV, each with unprece- 
dented precision Recently, the PAMELA collabora- 
tion has presented first results of the measurement of the 
positron fraction in the cosmic ray spectrum, which ap- 
pears to begin climbing quite rapidly between ~ 7 GeV 
and 100 GeV A similar trend was in fact also in- 
dicated by earlier experiments, including HEAT Q and 
AMS-01 [j, although with lesser statistical significance 
and over a smaller dynamical range. The behavior that 
seems to emerge in the positron fraction is very different 
from that predicted for secondary positrons produced in 
the collisions of cosmic ray nuclides with the inter-stellar 
medium (ISM). The situation is summarized in Fig [1] 
While unaccounted systematics in the measurements are 
in principle possible, we think it is worth reviewing what 
kind of physics may lead to such an energy spectrum; we 
shall argue that by far the simplest and most likely (as- 
tro)physical interpretation is that an additional, primary 
source of high energy positrons exists. 

The positron fraction is defined as 



f{E) 



(1) 



where the fluxes $i refer to the ones at the top of the 
atmosphere. Here and in the following, we shall keep 
implicit the dependence of fluxes on energy E. We re-cast 
the fluxes in terms of physically motivated contributions, 
without loss of generality: 

$e+ = V + S, (2) 

= A + V + S-V = ^^++A-V. (3) 

The term A represents the component of primary elec- 
trons accelerated in addition to any e"*" — e~ pairs; this 
term includes (but is not necessarily limited to) primary 
electrons accelerated in a typical ISM environment where 
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FIG. 1: The figure reports the positron fraction vs. energy 
measured by PAMELA 01, HEAT combined data 0] and 
AMS-01 A "typical" prediction based on purely secondary 
often used as a benchmark in the literature Q and power- 
law curves passing through PAMELA datum at 6.83 GeV 
with index e = 0.2 and e = 0.35 are also reported for illustra- 
tive purposes. 



no pairs are present. The term S represents the sec- 
ondary component of produced in hadronic cosmic 
ray collisions in the ISM. Note that an analogous term 
exists for e~ as well, which we denote S — V: there is 
indeed a small deficit V of secondary e~ compared to 
secondary e'^ due to charge asymmetry in the cosmic ray 
and ISM nuclei population, which is proton-dominated. 
Finally, we allow for a putative primary flux V of e^-e^ 
pairs, which is put to zero in typical predictions of f{E). 
A P-term might be due for example to unaccounted as- 
trophysical accelerators as pulsars or to a more exotic 
source as dark matter (DM) annihilation. In terms of 
these components. 



f{E) = 2 



A-V 

r + s 



(4) 



First, one trivial observation: since from particle physics 
we know that 2? > 0, from the empirical datum f{E) < 
1/2 it follows that A > V 0. This is not surpris- 
ing, since it is well known that accelerators of e~ (i.e. 
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electrons only, not pairs) exist in nature. Shocks on the 
background ISM at supernova remnants (SNRs, the most 
likely accelerators of Galactic Cosmic Rays) are the prime 
candidate in that sense; for modern predictions of the 
electron spectrum and an overview of past publications, 
see e.g. (6| and refs. therein. The low-energy behavior 
of /, typically found to decline from ~ 0.1 to ^ 0.05 up 
to 5 GeV, has been measured since longtime. This 
range is influenced by time and charge-dependent solar 
modulation, which is likely responsible for the differences 
among experiments. We shall neglect in the following so- 
lar effects which are irrelevant in the high energy range 
we focus on here. On the other hand, growing evidence 
has been collected in the recent years that f{E) might 
be rising at high energies, with the latest PAMELA data 
strongly favoring this observation. While one should wait 
for higher statistics and possibly an independent confir- 
mation (in particular by AMS-02 7]), it is useful to clas- 
sify the possible (astro)physical mechanisms leading to 
such an effect, a task we discuss in Sec. II. We will con- 
clude that the only one appearing viable requires V 0, 
which would imply the discovery of a new class of cosmic- 
ray sources (Sec. III). 



II. THE NECESSITY OF A PRIMARY 
SPECTRUM OF e+-e PAIRS 

A. Basic Arguments 

To prove the statement in the title of this Section, we 
shall adopt a "reductio ad absurdum" approach. Let us 
note that /' = df /dE > imphes 



A-V 
V + S 



< 0. 



(5) 



Let us now consider the Ansatz V — and neglect P for 
the moment (we shall see that this is justified, actually 
even a conservative assumption). Then, we should re- 
quire {A/Sy < to produce a rise. We see now that this 
requires highly implausible astrophysical conditions. 

In a very simple (position-independent) leaky-box 
model, the master equation for leptons simplifies into Q 

where $ is the lepton flux, q{E) the initial/injection spec- 
trum, t{E) oc E~^'' an effective containment time, and 
b{E) = —dE/dt ~ kE'^ the energy-loss rate function, 
which at the energies of interest (E^ 7 GeV) is dom- 
inated by synchrotron radiation and inverse Compton 
scattering. Let us assume spectra at the injection of 
the typical power-law form oc E^'^", E^'^f respectively 
for electrons and for protons (which, being the domi- 
nant hadronic component of cosmic rays and ISM, are the 
main responsible for secondary leptons). Protons suffer 
virtually no energy losses and obey an equation similar 



to Eq. ([6]) with b{E) — > 0, which yields for the steady 
state solution {d^/dt = 0) a spectrum oc E~'^p~^ with, 
a priori, 5 ^ 6e. The index 6 is constrained from the nu- 
clide ratio B/C to lie in the range 0.3 0.6 (a review of 
recent cosmic ray experiments and their interpretation is 
provided in ,9]). The convolution of this spectrum with 
the relevant cross-section is the source term for positrons, 
5+ oc E~^''~^ assuming an energy-independent inelastic 
cross-section and inelasticity. For the electrons, one has 
simply q_ oc E^'^" , plus a subdominant secondary con- 
tribution of similar magnitude and shape of the positron 
one. The resulting spectra of primary electrons and 
secondary positrons at the Earth are thus respectively 
oc E^"-, E-°'+ where a- = je+i, a+ = ■yp + i + S. Here 
£ symbolically represents the steepening due to diffusion 
and energy losses of leptons. For example, when one 
can neglect energy losses, $ oc q{E)T{E) and £ — Se', at 
sufhciently high-energy (TeV range) where energy-losses 
dominate it is easy to see that $ oc b{E)^^ J dE'q{E') 
and £ — \. Independently of the value of I, in this 
simple model we end-up with [A/S) oc _E~T'!+Tp+'' and 
thus Eq. (O requires 7p + (5 — 7e < 0. This condition 
seems extremely hard to achieve, requiring wildly differ- 
ent (by ^ 0.6) source spectral indexes for protons and 
electrons. This would contradict the standard theoreti- 
cal interpretation of the spectral difference between p and 
e~ observed at the Earth as due to different energy-loss 
properties, rather than intrinsic ones. Also note that the 
condition 7p -|- (5 — 7e < could not hold down to low 
energies, since the the sign of /' is negative around GeV 
energies. So, one should also invoke some spectral break 
in the injection electron spectrum placed ad hoc in the 
~ 7 GeV range. In summary, insisting in the prior V = Q 
and requiring thus {A/S)' < seems to imply: i) that our 
scenarios for the origin of Galactic electrons are wrong, 
requiring in turn cither new sources or new acceleration 
mechanisms different from the proton ones; ii) some de- 
gree of fine tuning, in the sense that the energy at which 
/' changes sign would correspond to some spectral break 
in the electron spectrum. These conditions appear way 
more extreme than allowing for primary sources of pairs, 
for which candidates (both astrophysical and exotic) do 
exist in the literature. 

On the top of the above considerations, there are em- 
pirical arguments which appear to disfavor this hypoth- 
esis. Defining ^tot = + ^e- i h'om Eq. (1) one has 



$tot(gl) ^ f{E2)^e+{El) 
<^tot{E2) f{Ei)<i>,^{E2) 



(7) 



Assuming only secondaries, the hardest spectrum the- 
oretically possible for positrons derived in jlQ] goes as 
~ £;-3.33 a^]-,ove 10 GeV. In the same paper it is reported 
that the softest possible spectrum for ^tot fitting (poorly, 
at 3 a) the data compiled in llj goes as ~ £;-3-54 ^^le 
same range. As a result, the maximal growth possible for 
the positron ratio — assuming only secondary positrons 
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and no priors for the electrons — is 



f{E2 



< 



E2 
El 



(8) 



As illustrated in Fig.[l] this is insufficient to fully explain 
the rise suggested by the PAMELA data. More in gen- 
eral, this argument proves that a better determination of 
the electron energy spectrum might reveal an inconsis- 
tency of the f{E) and ^tot data with a purely secondary 
origin of e"*", which does not resort to any theoretical 
considerations on the flux. 



B. Possible loopholes: further discussion 

Although illustrated in a simple leaky box scenario, the 
conclusion that 7-" is required appears robust. A sub- 
tle point in the considerations following Eq. (5) is that we 
really need to compare the spectrum of the primary elec- 
trons with the one of cosmic ray nuclei at much higher en- 
ergies (a factor ^ 20) since secondary leptons only carry 
a limited amount of the parent nucleus energy. Although 
a concavity of the spectrum would be naturally accom- 
modated in non- linear acceleration models (12j . there is 
no evidence that the (well measured) proton spectrum 
presents a noticeable change of slope around TeV energy, 
certainly not at the level of ~ 0.5. 

Another way around the previous conclusion may be 
to consider a progressively rising role of Helium nuclei 
as source of secondaries. Still, at energies around the 
TeV one should require a flux of Helium nuclei compa- 
rable to proton one and, at the same time, its spectral 
index harder than the proton one by an amount larger 
than 5. Actually, some indications of a hardening of He- 
lium spectrum has been claimed by the ATIC-2 collabo- 
ration [isj . But its amount and the energy range where 
it happens appear insufficient to explain the beaviour of 
f{E). For example, between 200 GeV and 1 TeV the flux 
ratio $p/$Hc varies by only ~ 15%. To exclude this pos- 
sibility, however, it would be important to compare the 
positron fraction and p and He spectra measured with 
the same instrument. Preliminary results by PAMELA, 
for example, do not support such an explanation since 
they show that both p and He fluxes are well fitted with 
the same spectral index ~ 2.73 up to 500 GeV jlj]. 

In 0, the possibility was discussed that the average 
interstellar proton spectrum may be harder than the one 
measured "locally" by an index ~ 0.15, invoking both a 
better agreement with the diffuse gamma-ray spectrum 
and the HEAT data on the positron fraction. Note that 
even this ad hoc adjustment would be insufficient to ex- 
plain the present evidence supporting /' > 0. Still, as- 
suming that this argument is correct, one would expect 
two qualitative predictions: since the "collecting cosmic 
ray volume" in the Galaxy depends both on primary type 
and energy (via diffusion and spallation effects) , a spatial 
non-universality of cosmic ray acceleration should reflect 



into a change of slope of a single species vs. energy, and 
of different species from one another. The TRACER col- 
laboration has instead reported that for nuclei between 
Oxygen and Iron a single power-law index~ 2.7 can fit 
all the data in the GeV to TeV energy per nucleon, with 
possible variations within ~ 0.05 [la|. We take here the 
agnostic point of view that the gamma-ray "excess" is 
not understood at the moment (it might even be due to 
a calibration problem, see [HI), but note that it may be 
considered as well as an indication that additional sources 
exist contributing above the GeV range, which is consis- 
tent with the hypothesis of additional primary emitters 
of cosmic ray positrons. 

One may further wonder if a rising positron ratio might 
be due to an unexpected energy-dependent beaviour of 
the diffusion index; from previous considerations and in 
the simplest case of 7p = 7e, it would follow indeed 
{A/S) oc . For the sake of the argument and with 
a slight abuse, let us assume 5 to be "slightly" energy 
dependent (this is not rigorous since the previous solu- 
tion has been derived for a constant 5.) The condition 
{A/Sy < translates into the requirement that S{E) 
declines with energy faster than l/\n{E) in the 0.1-1 
TeV energy range of the parent nucleus producing the 
positrons. This is a relatively large effect, with 6 drop- 
ping by at least a factor 3 in a decade of energy to 
account for the rising f{E). Even the proposal that i5 
changes from ~ 0.6 above 10 GeV to ^ 0.3 at TeV ener- 
gies ( [6] and rcfs. therein) appears insufficient to account 
for the sign of /'. On the other hand, this argument 
faces another difficulty: the featureless power-law of CR 
protons would result from a fine-tuned compensation of 
the variation of S{E) and the injection spectrum, which 
seems improbable, the two being unrelated. It is worth 
noting however that even this baroque scenario is testable 
empirically from high-energy B/C data. 

Another approximation in the previous argument is the 
assumption of a constant cross-section. A rising inelas- 
tic cross section would reflect in the secondary energy 
spectrum. Indeed, the inelastic cross section grows with 
energy (see e.g. [i3l), but only logarithmically, by ~ 1% 
between 10 and 100 GeV and by ~ 10% between 0.1 and 
1 TeV, i.e. equivalent at most to a power-law of index 
^ 0.04. This is more than one order of magnitude smaller 
than what needed to explain the positron feature. 

Finally, let us come back to D, or better "D/S. The 
above considerations are a fortiori true if {'D/S)' < 0, i.e. 
if the relative difference between positrons and electrons 
remains constant or declines with energy. Note that 
this function is mostly dependent on particle physics, 
apart for the convolution with the primary cosmic ray 
spectrum. If wc take for example the e'^ and yields 
for a proton power-law spectrum reported in Fig. 12 of 
Ref. [13] one can conclude that: i) the secondary spectra 
are slightly harder than the primary one, consistently 
with the energy-dependence effect discussed above; 
ii) the electron spectrum is slightly harder than the 
positron one, i.e. V/S is slightly decreasing with energy. 
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Although a detailed study would be required to assess 
more quantitatively the uncertainties in this argument, 
it is clear that a sufficiently strong dependence of V/S 
which might account for a rise in f{E) appears out of 
question. 



III. CONCLUSIONS 

In summary, motivated by observational evidence and 
in particular by PAMELA data [3] , we have discussed un- 
der which conditions a rise in the positron fraction f{E) 
can take place. Barring the case of systematics in the 
measurements, we have analyzed the following hypothe- 
ses: 

1) "Anomalous" primary electron source spectrum. 

2) Spectral feature in the proton flux responsible for 
the secondaries. 

3) Role of Helium nuclei in secondary production. 

4) Difference between local and ISM spectrum of pro- 
tons. 

5) "Anomalous" energy-dependent behaviour of the 
diffusion coefficient. 

6) Rising cross section at high energy. 

7) High energy behavior of the e^/e~ ratio of secon- 
daries h\ pp collisions. 

All of the above options seem to be at least strongly 
disfavored if not already ruled out; nonetheless, we have 



summarized the signatures associated with different 
explanations, and the way to test them observationally. 
Among the options listed which assume V — the 
one coming closer to a (very bad) fit to the data is 
nr. 1), which is not only disfavored by the data, but 
requires an ad hoc adjustment, lacking at the moment 
an astrophysical model producing it. We concluded 
that the most likely cause of the energy trend of the 
positron fraction is the presence of a primary flux of e^, 
which — both in astrophysical and exotic models — are 
probably injected in the form of pairs (see for an 
early review of possible primary sources) . Accepting this 
solution, one has at high energies f{E) ~ (2 -I- ^)~^; 
from a rise at > 7 GeV one can further deduce that 
the spectrum of pairs is harder than the one of primary 
electrons, which is also a typical prediction in pulsar or 
DM annihilation/decay models. In the opinion of the 
author, the positron spectral shape and normalization 
suggest pulsars as the most plausible responsible for the 
emission, a possibility which has drawn some attention 
lately [1^, [20, [21| . At very least, these objects should be 
seriously considered among the main actors of the high 
energy Galactic sky; perhaps they are also responsible 
for most unidentifled Galactic gamma ray sources, as 
recently supported by the discovery by FERMI of a 
pulsating gamma-ray emission from the SNR CTA 1 . 
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